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ABSTRACT 

We present an analysis of approximately 13-yr of observations of the intermittent pulsar 
B 193 1+24 to further elucidate its behaviour. We find that while the source exhibits a wide 
range of nulling (~ 4 — 39 d) and radio-emitting (~ 1 — 19 d) timescales, it cycles be- 
tween its different emission phases over an average timescale of approximately 38 d, which 
is remarkably stable over many years. On average, the neutron star is found to be radio- 
emitting for 26 ± 6 % of the time. No evidence is obtained to suggest that the pulsar under- 
goes any systematic, intrinsic variations in pulse intensity during the radio-emitting phases. 
In addition, we find no evidence for any correlation between the length of consecutive emis- 
sion phases. An analysis of the rotational behaviour of the source shows that it consistently 
assumes the same spin-down rates, i.e. v = —16 ± 1 x 10~ 15 s~ 2 when emitting and 
v = —10.8 ± 0.4 x 10~ 15 s~ 2 when not emitting, over the entire observation span. Cou- 
pled with the stable switching timescale, this implies that the pulsar retains a high degree 
of magnetospheric memory, and stability, in spite of comparatively rapid (<~ ms) dynamical 
plasma timescales. While this provides further evidence to suggest that the behaviour of the 
neutron star is governed by magnetospheric-state switching, the underlying trigger mechanism 
remains illusive. This should be elucidated by future surveys with next generation telescopes 
such as LOFAR, MeerKAT and the SKA, which should detect similar sources and provide 
more clues to how their radio emission is regulated. 

Key words: pulsars: individual: PSR B1931+24 - pulsars: general. 



1 INTRODUCTION 

PSR B 193 1+24 was di scovered in 1985 w ith the NRAO 100-m 
Green Bank Telescope dstokes et al.lll985h . However, it was not 
until 13 years later, during routine pulsar timing observations at 
Jodrell Bank, that the remarkable properties of this object were 
noted. Unlike conventional nulling pulsars, which undergo tem- 
porary emission cessation (or pulse nu l ling) over timescales of 
< 100 pulse periods (e .g. lBackerlll970l lRitchings|[l97o1 : iRankml 
1 19861 : IWang et alj|2007b . PSR B 193 1+24 has been found to ex- 
hibit extremely long duration, active (~ 5 — 10 d) and quiescent 
(~ 25 — 35 d) radio emission phases. Remarkably, these emis- 
sion phases also repeat quasi-periodically, are broadband in fre- 
quency and are found to be correlated with the rotational behaviour 
of the star. That is, the source exhibits a spin-down rate which is 
~ 50 % greater in the active (radio-on, hereafter) emission phases, 
compar ed with that of the n on-radio emitting (radio-off, hereafter) 
phases dKramer et al . 2006). 
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To date, a handfu l of similar objects have been discovered, 
e.g. PSR 1183 2 +002 9 dLorimer et ai]|2012h and PSR J1841-0500 
dCamilo et al .1 |2012|) . Despite this, only the properties of 
PSR B 193 1+24 have been studied, or attempted to be explained, in 
some detail (see. e.g.jKramer et a l. 2006; C ordes & Shannonl2008l : 
iRea et alj2008l ; lRosen et alj201 lklLi et al.l2012l) . As such, theories 
which attempt to explain the behaviour of such long-term intermit- 
tent pulsars are centred mainly on observations of this object, par- 
ticularly those in the radio regim^B 

Amongst the most promising of these theories, is the pos- 
sibility that the pulsar magnetosphere experiences systematic 
variation in its global charge distribution; that is, it under- 
goes 'magnetosphe r ic-state switching ' (e.g. B artel et all 19821; 
Contopoulosl 120051: iLvne et all l201ol ; iTimokhinl l20ld : Li et al .1 



2012b . 



In this context, it is proposed that global re-distributions 



1 While PSR B 1 93 1 +24 has been obs erved both in the radio IStokes et all 
ll985tlKramer et al.ll2003) and optical <Rea et alj|20O8b regimes, it has so 
far only been detected in radio observations. 
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of current are responsible for causing alterations to the morphol- 
ogy of the emitting region, which result in the presence (or ab- 
sence) of radio emission, as well as fluctuations in the spin- 
down rate via associated changes in the magnetic field config- 
uration. However, it is unclear how these alterations occur, nor 
what their intrinsic timescales should be; an array of vastly dif- 
ferent tr igger mechanisms have been proposed, e .g. orbital com- 
panions l lRea et al . 2008; Cordes & Shannon 2008]), non-radial os - 
cillations dRosen et alj|201 if), precessional torques {jonesl 1201 lb - 
magnetic field instabili ties dGeppert et al .120031 ; lUrpin & Gill2004l ; 
iRhemhardt et alj|2004) and polar cap surface temperature varia- 
tions dZhang et alj|l99% . but no consensus has been reached on a 
definitive mechanism. Subsequently, little is known about the pro- 
cesses which govern the behaviour of intermittent pulsars, nor what 
their relationship is with conventional nulling pu lsars and rotating 
radio transients (RRATs : lMcLaughlin et alj|2006h . 

In this paper, we use an unparallelled span of observations 
to probe the emission and rotational properties of PS R B1931+24 
in deta il. We focus on expanding upon the results of I Kramer et al.l 
j2006h . by using our longer observation baseline to investigate 
both short- and long-term variations in the behaviour of the pulsar. 
We present the details of the radio observations in the following 
section. In Section [3] we describe the emission variability of the 
source. This is followed by an investigation of its rotational stabil- 
ity in Section [4] We discuss the main implications of our findings 
in Section[5]and present our conclusions in Section[6] 



2 OBSERVATIONS 

In this work, we present an approximately 13-yr span of observa- 
tions (29 April 1998 to 19 May 20 11), including previously pub- 
lished data (see lKramer et alj2006l) . which has been used to inves- 
tigate the emission and rotational properties of PSR B1931+24. 
These data were predominantly obtained with the 76-m Lovell 
Telescope and the 28x25-m Mark II Telescope at Jodrell Bank. 
Data was also obtained with the 94-m Nancay Radio Telescope 
in France, so as to bridge a gap in observations during a period 
of extended telescope maintenancfl Two back-ends were used to 
acquire the Lovell observations: the Analogue Filter Bank (AFB; 
up to May 2010) and the Digital Filter Bank (DFB; since Jan- 
uary 2009). Table[T]shows the typical observing characteristics for 
these instruments and the Nancay data. We note that while the av- 
erage observing cadence is less than once per day for the AFB, a 
series of more intense observations, i.e. approximately twice-daily 
monitoring, was initiated from 2006 onwards to provide better con- 
straints on the emission phase transition times. 



3 EMISSION VARIABILITY 
3.1 Pulse intensity modulation 

Following the work of iKramer et al. I d2006f) . we sought to eluci- 
date the emission variability in PSR B 193 1+24, based on a longer 
observation baseline and higher cadence data. In particular, we 

2 Between 20 October 2004 and 2 February 2005, the Lovell Telescope was 
unable to observe PSR B1931+24 due to the discovery and replacement 
of a cracked tire. The Mark II Telescope was used for continued timing 
measurements of other sources during this time. 



Table 1. System characteristics for observations of PSR B1931+24 from 
29 April 1998 to 19 May 201 1 . The total time span and number of observa- 
tions are denoted by T and N respectively. The most common observation 
duration is represented by At and the typical cadence is denoted by C. The 
typical centre frequency, bandwidth and channel bandwidth of the observa- 
tions are given by u, B and B c ^ respectively. 



System Property 


AFB 


DFB 


Mark II 


Nancay 


T(d) 


4403 


856 


3661 


243 


N 


3996 


1168 


636 


78 


At (min) 


12 


12 


42 


13 


f(d 1 I 


0.9 


1.4 


0.2 


0.3 


i^(MHz) 


1402 


1520 


1396 


1368 


B (MHz) 


32 


384 


32 


64 


B ch (MHz) .... 


1 


0.5 


1 


4 



were interested in further characterising the timescales of emis- 
sion variation, and the durations of the radio-on and -off emission 
phases. For the purpose of t his work, we used the activity duty cycle 
(ADC) method discussed in lYoung et all d2012h . Here, we visually 
inspected the average pulse profile data, which was formed over the 
entire integration length and bandwidth of observations, and con- 
structed a time-series of one-bit data corresponding to the 'radio 
activity' of the pulsar; that is, 1 's for profiles consistent with the 
emission signature of the source and 0's for noise dominated pro- 
files (i.e. radio-off states). The result of this analysis for the highest 
cadence observations (i.e. the last ~ 4 yr of data) is shown in Fig.Q] 
In the data set, some of the observations which define the 
timescales of emission are separated by several days, particularly at 
earlier epochs. Therefore, a number of the inferred emission phase 
durations in the data do not accurately represent the behaviour of 
the pulsar. Consequently, we only consider emission phase dura- 
tions which have gaps < 5 d between consecutive observations, so 
that we can pinpoint transition times (taken to be the midpoints be- 
tween radio-on and -off observations) as accurately as possible. Af- 
ter applying these criteria we were left with only the highest confi- 
dence emission phase durations, which are shown in Fig. [2] We note 
that there is no evidence for emission cessation over timescales less 
than a day. The average times in each emission phase are 8 ± 4 d 
and 22 ± 7 d for the radio-on and -off phases respectively. These 
values are c onsistent with the av erage radio-on and -off timescales 
observed in IKramer et all d2006h . that is ~ 6 d and ~ 28 d respec- 
tively. We note here that the radio-on and -off timescales have a 
wi der range of values than previously thought (Fig. [2] c.f. Fig. lc 
of IKramer e"tai]l2006l) . The improvement in the determination of 
these values is due to the increased data span and observation ca- 
dence. 

We used these same emission phase durations to calculate the 
ADC; that is, the percentage of time spent in the radio-on phase, 
compared with the total observation time. Given that there are gaps 
between consecutive emission phases, and that there is a disparity 
between the number of suitable radio-on and -off phases, we deter- 
mined the average total time spent in each emission phase through 
bootstrapping the observed distributions of the emission phase du- 
rations. That is, the data was sampled with replacement to obtain 
10 6 resamples, and average durations, for each distribution of emis- 
sion phases. The resultant values were used to calculate the average 
ADC over the entire data-set. The error in the ADC was calculated 
from the uncertainty in the total time spent in the radio-on phase, 
which was obtained from bootstrapping the observed uncertainties. 
This results in an ADC = 26 ± 6 % which is larger than, but 
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Figure 1. Sequence of observations, denoted by the black lines, carried out over ~ 4 yr (12 January 2007 to 19 May 2011). The data is separated into 5 
contiguous 318-d panels. The times of observation and the times when PSR B 193 1+24 was radio-on (full height) and radio-off (half height) are shown by the 
extent of these lines. 
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Figure 2. Histograms showing the durations when PSR B 193 1+24 is observed to be in a particular emission phase, radio-on (left) and radio-off (right). 



consistent with, the ItCramer et al.l d2006h value of 19 + 5%. We 
note that the larger uncertainty in our value most likely results from 
slight changes in observation cadence throughout the ~ 13-yr data 
set. 

As our average ADC is the same as in l Kramer et al. I fcood) and 
yet the range of radio-on and -off durations is wider, we wanted to 
ascertain if there was any temporal ev olution in these p roperties. 
Following the stride fitting method of iLvne et alj ( l2010h . we de- 



termined the ADC over segments of length T = 100 cjf], offset 
by intervals of T/4 = 25 d across the data set. We initially esti- 
mate the uncertainty in ADC for each data segment using the sum 
of the errors in the transition times between consecutive emission 

3 The interval length T was chosen to provide a compromise between res- 
olution in the ADC and our sensitivity to short-term noise variations, that 
is short-term fluctuations which do not reflect the typical behaviour of the 
object. 
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Figure 3. The evolution of the activity duty cycle of PSR B1931+24 over 
approximately 13 yr of data. 



phase durations. Due to the irregular time sampling, however, there 
are also gaps between consecutive observations during emission 
phases. As the minimum time spent in a consecutive radio-on and 
-off phase is approximately 10 d, we assume that any gap between a 
neighbouring observation which exceeds this duration should con- 
tribute to the total uncertainty in the ADC. In order to compromise 
between data accuracy and volume, we apply cut-offs to the un- 
certainties in the data; that is, we do not consider intervals which 
contain observations separated by more than 25 d or whose error in 
the ADC is greater than 50 %. The resultant ADC as a function of 
time is shown in Fig. [3] 

We average these values (Fig. |3j to obtain 
(ADC) = 26 ± 7%, which is consistent with that obtained 
by bootstrapping the emission durations. We attribute the signifi- 
cant error in this quantity to the notable fractional uncertainties in 
the ADC throu ghout the data set . By implementing an Anderson- 
Darling tes{3 !; PressetalJll992h we find that the ADC values 
are consistent with coming from a normal distribution, assuming 
Gaussian statistics for the uncertainties. There is, therefore, no 
compelling evidence for significant variation in the ADC of 
PSR B 193 1+24. 

We also sought to determine whether there were any interde- 
pendencies in the characteristics of the pulsar emission. In partic- 
ular, we were interested in discovering whether there is a correla- 
tion between the length of successive radio-on and -off phases. For 
this analysis, we selected nine data intervals from the entire data 
set which had the highest observation cadence and, thus, the high- 
est confidence in emission phase durations. The properties of these 
data intervals are listed in Table [2] We correlated the durations of 
successive radio-on and -off emission phases, for each of these data 
intervals, to determine whether there was any connection between 
the length of consecutive emission phases (see Fig.|4j- 

We find that there is no significant correlation between the 
length of time in a given emission phase and that of the opposing 
mode consecutive to it. This suggests that the pulsar does not retain 
a memory of the length of its previous magnetospheric-state after 
a transition. We also computed the ADC for each data interval to 
determine whether the results of the previous analysis may have 



4 This test is adapted from the Kolmogoro v-Smimov test to have greater 
sensitivity towards the tails of a distribution ( Press et al. 1992), which make 
it better suited to testing for normality. 




Figure 4. The correlation between the radio-on and -off emission phase du- 
rations of PSR B 1 93 1 +24. Emission duration histograms for all the radio- 
on (top panel) and radio-off (middle panel) phases from the observation 
periods defined in Table [2] as well as the linear regression of the two pa- 
rameters (bottom panel). There is no evidence for a correlation between 
consecutive radio-on and radio-off interval lengths. 



been biased by observation cadence (see Table [2j. We again find 
that there is no significant evidence for variation in the ADC over 
time. 

We were also interested in determining the average flux den- 
sity of the pulsar, attributed to the dif ferent phases of emission , 
using the modified radiometer equation dLorimer & Kramej2005h : 



S = 



l3SNRT sys 



G^BT]j P 



(1) 



Here, /3 ~ 1 is the digitisation factor, G ~ 1 Jy K _1 is the tele- 
scope gain, T S ys ~ 35 K is the average system temperature, n p — 2 
is the number of polarisations, B — 384 MHz is the observing 
bandwidth, P — 814 ms is the pulsar period and W cq = 14 ms is 
the equivalent pulse width. In total, we averaged 211 radio-on and 
106 radio-off DFB observations (see Fig. [5}, to obtain total integra- 
tion times (T) of approximately 40 hr and 20 hr respectively. Con- 
sequently, we place a limit on the mean flux density in the radio-off 
phase 5*off < 2.0+0.4 /ijy (assuming a limiting SNR ~ 3). For the 
radio-on phase, we obtain a mean flux density S on =40 + 8 /xJy 
(SNR ~ 90), which is at least 20 times brighter than that in the 
radio-off phase. 

Converting th e flux densities into pseud o-luminosities, using 
Z/1400 = 51400 d dLorimer & Krameill2005h . we find L1400, off < 
0.04 mjy kpc 2 and L1400, on ~ 0.84 mjy kpc 2 ( for a pulsar dis- 
tance d ~ 4.6 kpc, from the NE2001 model; ICordes & Laziol 
120021) . We note that the upper limit on the pseudo-luminosity of 
PSR B 193 1+24 in the radio-off phase is barely brighter than the 
two weakest kn own radio pulsars i.e . PSR J0030+0451 (£1400 ~ 



0.04 mjy kpc 2 ; iLommen et alJl200Ch and J2144-3933 (Li.„. 
0.02 mjy kpc 2 ; iLorimeil Il994 ). This implies that the radio-off 



phases are consistent with emission cessation. However, it may 
be that the source exhibits extremely weak, underlying emission 
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Table 2. Properties of the nine data intervals chosen to characterise the correlation between the radio-on and -off emission phase durations of PSR B 193 1+24. 
The start and finish times of each data interval, of length T, are denoted by MJD s tart and MJDfj n ; s h respectively. The number of emission phases which are 
radio-on and -off are given by N on and iV Q g respectively. The average total time spent in a radio-on and -off phase are t on and t g , and the activity duty cycle 
of the pulsar for each data interval is denoted by ADC. 



MJD atart 


MJD finiah 


T(d) 


AT on 


N oB 


*on (d) 


t oft (d) 


ADC (%) 


51814.07 


51924.38 


110.31 


4 


3 


7±1 


24 ±2 


22 ±8 


52762.52 


52921.53 


159.01 


5 


4 


10 ±4 


27 ±4 


28 ±2 


53733.23 


53832.01 


98.78 


4 


3 


7± 2 


22 ±7 


29 ±5 


53889.35 


54084.26 


194.91 


8 


7 


5± 1 


21 ±3 


19 ±5 


54176.30 


54504.73 


328.43 


11 


10 


9± 1 


23 ±2 


27 ±5 


54609.06 


54743.66 


134.60 


5 


4 


8±1 


23 ±1 


26 ±4 


54775.96 


54973.34 


197.38 


6 


5 


11 ± 2 


26 ±4 


30 ±4 


55066.68 


55489.50 


422.82 


12 


12 


7± 1 


28 ±2 


21 ±3 


55563.29 


55688.01 


124.72 


4 


4 


12 ±1 


19 ±6 


38 ±3 
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Figure 5. Average pulse profiles of PSR B 193 1+24 showing the charac- 
teristics of the radio-on (top) and radio-off emission (bottom). The peak 
pulse intensities are normalised to one for each emission mode profile. The 
radio-on profile is offset in intensity from the radio-off profile for clarity. 
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Figure 6. The evolution in the ratio of the first peak to second peak inten- 
sities of PSR B 193 1+24, over ~ 165 d of DFB observations. The average 
and ±lcr values are denoted by the (dashed) and (dot-dashed) lines re- 
spectively. The fluctuation in the peak ratio values does not appear to be 
significant with respect to the average. 



which is far below our detection threshold. This could be com- 
pared to PSR B0823+26, another intermittent radio source, which 
exhibits a factor of 100 difference i n in tensity between its se parate 
emission modes dYoung etalj|2012l c.f. lEsamdin et alj|2005t) . 

Observations of PSR B1931+24 bv lKramer etall ( feOOfj) sug- 
gest that it exhibits enhanced particle flow during its radio-on 
phases. We were interested in seeing whether other, less severe, 
changes in particle flow occurred which might also be attributed 
to the mechanism which determines its emission behaviour. One 
way to test whether the emission in these phases was consistent 
with a 'steady- state ' particle flow (c.f. Crab enhanced emission, see 
IShearer et al]2003l) . was to determine whether the pulsar exhibited 
any systematic pulse intensity fluctuations during radio-on phases. 
Such analysis is ideally performed using single-pulse data which, 
in this instance, was not available. Therefore, we used the 12-min 
average pulse profiles to determine whether there were any vari- 
ations in pulse intensity, which might correlate with the position 
of an observation around an emission phase transition (c.f. pulse 
intensity decay and rise times around nulls in PS R B0809+74; 
iLvne & Ashworthl [l983l ; Ivan Leeuwen et al.ll2002h . We find that 
there is no significant correlation between the pulse intensity prior 
to, or after, a radio-off phase. We also find that there is no signifi- 
cant modulation in the emission during a given radio-on phase. That 
is, the variations in pulse intensity are dominated by random fluc- 



tuations. This suggests that the particle flow in the magnetosphere 
of PSR B 193 1+24 remains constant during the radio-on emission 
phases, to the limit of our measurement sensitivity and the intrinsic 
flux variation of the source. 

To check whether there were any pulse shape changes dur- 
ing the radio-on ph ases, as seen in ot her pulsars exhibiting period 
derivative changes l lLvneetal.ll2010l) . we examined the variation 
in the ratio of the first to second component peak intensities over 
time. For this analysis, we aligned DFB profiles with an analytic 
templat^fl using a \ 2 minimisation method to match the pulse lon- 
gitudes of the peak bins. This, in turn, enabled us to quantify the 
component peak intensities over time. To obtain a correct represen- 
tation of the peak-to-peak ratios, we only performed this analysis 
on the highest SNR observations (SNR > 10). Figure [6] shows the 
distribution of peak ratios. The average peak ratio is 0.69 ± 0.09 
and is found to be consistent with that of the analytic template (i.e. 
~ 0.7). 

In order to determine whether the peak ratio variations were 
significant, we again performed an Anderson-Darling test on the 



5 This template was produced with the paas program, which was used to 
fit von-Mises functions to the highest SNR observation. For an overview 
see http://psrchive.sourceforge.net/changes/v5.0/. 
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distribution of values, and found no evidence to suggest the data 
departs from a normal distribution, assuming Gaussian error bars. 
As a consistency check, we also implemented a reduced \ 2 test to 
quantify the fluctuations in the pulse shape. For these observations 
we obtain \ 2 ~ 0.92, which is consistent with random noise dom- 
inating the profile variation. 

3.2 Periodicity analysis 

To co mplement, and improve upon, the results of iKramer et al.l 
d2006h we have conducted an extensive analysis on the periodic- 
ity of radio emission from PSR B 193 1+24, using a much longer 
data set (i.e. ~ 13 yr of data). Here, we use a modified wavelet 
analysis, i.e. weighted w avelet Z-statistic (WWZ) analysis dFosted 
ll996UYoungetalJl2012l) . to reveal any periodic variations within 
the data. 

For the purpose of our analysis, we chose a wavelet tuning 
constant c = 0.001 so that we could strike a balance between 
frequency and time resolution. The resultant WWZ transform of 
the one-bit time-series data, showing the evolution of the spectral 
power at successive epochs (or time lags), is displayed in Fig. [7] In 
the 2-D transform plot, we note that the peak frequencies, i.e. WWZ 
frequency maxima, modulate over time. This, in turn, gives rise to 
several features in the integrated spectrum. The broad prominent 
peak in the integrated spectrum shows that the WWZ power is pref- 
erentially distributed around ~ 0.031 - 0.02 4 d" 1 (~ 32 - 42 d) , 
which corresponds well with the results of IKramer et alj J2006h . 
We also note the presence of other features in the transform plane 
at ~ 0.021 d _1 (~ 48 d), ~ 0.035 d _1 (~ 29 d), approximately 
0.042 d" 1 (~ 24 d) and ~ 0.052 d" 1 (~ 19 d) over small time 
intervals. However, they generally have much lower spectral power 
which casts doubt on their significance. We find that the features 
centred on ~ 0.042 d _1 and ~ 0.052 d _1 are most likely spuri- 
ous due to their proximity to a long, approximately 100-d gap in 
the datfl We note that the better observation sampling towards 
later times (after MJD ~ 52810) results in better resolution of 
the intrinsic variation of the source, that is the fluctuation peri- 
ods are generally represented by greater WWZ power towards later 
tim es. We also conside r the error estimation of these data, follow- 
ing lYoungetalJj2012l) . in the next subsections. 

3.2.1 Confusion limit estimation method 

To estimate the significance of spectral f eatures, we first imple- 
mented the confu sion limit method (see iTempleton et alj 1 20051 ; 
lYoung etal.ll2012l) . For this analysis, we located the peak fluc- 
tuation frequency (taken to be the frequency bin with the great- 
est spectral power) at each epoch and estimated its maximum 1-a 
uncertainty as the half-width at half-maximum of the Z-statistic, 
Z(oj, t). The result of this analysis is shown in Fig. [8] We note 
that the peak period spans ~ 20 — 50 d over time, with a few 
epochs showing periodicities of ~ 100 d. However, the features 
in the wavelet transform attributed to this longer period are only 
represented by low WWZ power, and are over an interval of time 
which has a number of gaps between successive epochs (MJD ~ 
51987 - 52034, ~ 52040 - 52077 and ~ 52207 - 52240). There- 
fore, we do not believe that the longer fluctuation period is intrinsic 
to the source. 

6 Nothing meaningful can be obtained for a fluctuation period which is 
shorter than the length of a gap between data points. 
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Figure 8. Peak WWZ fluctuation frequencies (crosses) and periods (open 
circles) for the ~ 13-yr PSR B1931+24 radio emission activity data. Error 
bars are l-cr values computed using the confusion limit estimation method. 
The average error in the fluctuation frequency (or period) is approximately 
6 %. Note that the peak fluctuation period spans ~ 20 — 50 d apart from 
at a few epochs where the data sampling is likely to have affected the local 
matching. 



We performed another Anderson-Darling test on the peak- 
fluctuation period data to determine whether the variations were 
significant. From the results of this test, we find that the short-term 
modulation in the source's periodicity (see steps in Fig. [8} is not 
consistent with random variation. In fact, we find that this modu- 
lation becomes more significant with increasing observation sam- 
pling. As a result, we attribute these short-term fluctuations to the 
quasi-periodic nature of the object. 



3.2.2 Data-windowing method 

As a consistency check, we also used a data-windowing method to 
compute a measure of the m odulation in the peak frequency and pe- 
riod (see lYoung et al.l2012h . Here, we split the WWZ data into 400- 
d segments and, for each of these segments, determined the me- 
dian fluctuation frequency and corresponding period (see Table[3]l. 
These quantities were then compared with the 'noise background' 
of the WWZ data to determine their significance. Here, the noise 
background was estimated by computing the standard deviation of 
the transform data, by bootstrapping over several hundred realisa- 
tions. Subsequently, we find that the first four data segments do not 
have median peak WWZ values greater than the 5a ~ 60 signif- 
icance level derived from bootstrapping. However, the rest of the 
segments (excluding the data spanning MJD = 54150 — 54550) 
have values greater or equal to this cut-off, which indicate that 
the periodicities within these data intervals are intrinsic to the 
source. We note that the uncertainty in the peak fluctuation fre- 
quency becomes consistently smaller at later epochs, which results 
from an increase in resolution due to the improved data sampling. 
From this analysis, we obtain median peak fluctuation frequencies 
~ 0.035 - 0.021 d _1 (~ 28 - 47 d), which are largely consistent 
with the previous analysis. Through performing another Anderson- 
Darling test, we find that the variations in the long-term periodicity 
of the source, i.e. over several years, are not significant. As a re- 
sult, we infer that the overall average periodicity in the pulsar radio 
emission (38±5 d) is highly stable over long timescales (i.e. years). 
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Figure 7. WWZ transform of the ~ 13-yr PSR B 193 1+24 radio emission activity data-set (left) and the corresponding integrated power spectrum (right). The 
integrated power spectrum reaches a maximum at ~ 0.028 d _1 (~ 36 d). Throughout the transform plane the peak frequency modulates over time, typically 
ranging ~ 0.024 - 0.032 d~ x (~ 31 - 42 d). The bootstrap 5-ct significance level con'esponds to Z(ui, r) fa 60 for this data-set. 



Table 3. Results from the data-windowing WWZ error analysis. We quote 
the MJD ranges of the data windows analysed, the corresponding median- 
peak WWZ values WWZ max, as well as the peak fluctuation frequencies 
ffl uc and periods Pfl uc . The standard 1-cr uncertainties are quoted in the 
parentheses and are in units of the least significant digit. 



MJD range 


wwz max 


u auc (W- 2 d" 1 ) 


Pfluc(d) 


50950-51350 


18(3) 


3.2(6) 


31(6) 


51350-51750 


30 (10) 


2.5 (3) 


40(5) 


51750-52150 


28 (7) 


3(1) 


30 (40) 


52150-52550 


23(4) 


2.6 (7) 


40 (30) 


52550-52950 


60 (20) 


2.1 (4) 


47(7) 


52950-53350 


60 (20) 


3.0(7) 


34(6) 


53350-53750 


90 (10) 


2.75 (4) 


36(1) 


53750-54150 


80 (30) 


2.7 (3) 


37 (3) 


54150-54550 


50 (20) 


3.5 (5) 


28 (4) 


54550-54950 


80 (30) 


2.80 (5) 


36(1) 


54950-55350 


70 (30) 


3.1(2) 


32(3) 


55350-55700 


100 (20) 


2.4205 (1) 


41.3 (2) 



3.2.3 WWZ analysis summary 

The results from the WWZ analysis above show that 
PSR B 193 1+24 exhibits quasi-periodicity in its radio emis- 
sion switching over timescales of weeks to months. That is, 
the source displays a relatively wide range of periodicities 
(~ 20 — 50 d), when considering its short-term emission variation. 
Over longer timescales (i.e. years), however, the neutron star 
exhibits a highly stable average periodicity (38 ± 5 d) in its radio 
emission. This implies that the mechanism governing the source's 



behaviour is a non-random systematic process, which is perturbed 
over timescales of weeks to months. 



4 SPIN EVOLUTION 

4.1 Residual fitting and charge density estimation 

To gain context on the emission modulation of PSR B 193 1+24, 
we studied its rotational behaviour. We were particularly interested 
in determining whether the object experiences any temporal evolu- 
tion in the spin-down rates associated with the different phases of 
emission which, in turn, may be correlated with alterations to the 
charge density d istribution of the pulsar magnetosphere (see, e.g., 
lTimokhij|2Qlch . 

To quantify the spin-down rate variation in PSR B 193 1+24, 
we took a different approach comp ared with convent ional phase- 
coherent timing analysis (see, e.g.. iLvne et alj|l996j). d ue to the 
emission cessation in this source. Kra mer et al. I J2006T) demon- 
strated the success of using a dual spin-down rate, least-squares 
fitting method to derive the spin-parameters for this pulsar over a 
short span of data and, hence, the spin-down rates corresponding to 
each emission phase. We have expanded upon this initial analysis 
by applying the same fitting proced ure to eight well-sam pled data 
sets (including the one presented in lKramer et al.|[2006h - The pro- 
cedure entails minimising the timing residuals of PSR B 193 1+24 
(calculated using a fixed i> ephemeris) by fitting for the change in 
rotational phase, period and period derivative (A<fl AP and AP 



7 The change in rotational phase is the offset required, with respect to an 
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Figure 9. Least-squares, weighted fit to timing residuals of PSR B 1 93 1 +24 
for observations ranging MJD ~ 54777 — 54973. The timing residuals 
(filled circles) here are calculated with respect to a fixed i> ephemeris, while 
the fitted model (solid line) assumes two spin-down rates (see Eqn. Ir- 



respectively) with respect to the separate emission phases. The val- 
ues of P onj g, and hence z> on> g, were derived from the difference 
between the average spin-down rate (P av ) and AP associated with 
these modes. The total A<j>, AP and AP, which were required to 
obtain a phase-coherent solution, were modelled via: 



tr 



= A,.P + ^ 



(t-ti) 



AP . 
+ 2P •(*■ 



(2) 



where t res and t\ are the residual (t bs — iprcd) and reference times 
respectively (here £ pr cd represents the time-of-arrival predicted by 
a single spin-down rate model), and AP; is the change in spin-down 
rate in each emission phase (i.e. ±AP). An example of this fitting 
process is shown in Fig. [9] The rotational behaviour of the pulsar is 
clearly characterised well by the model. 

We note that the total contribution of !> on , g to the aver- 
age is determined from the total time spent in a given emission 
phase, which is inherently dependent on the transition times be- 
tween phases. Typically, we can only constrain transition times 
to within about a day which, in turn, imparts systematic errors 
into this analysis. To propagate the effects of such errors we 
performed a Monte-Carlo simulation on the data. Here we per- 
formed 10 fits, for each data interval, using transition times 
which are randomly distributed between the last radio-on and first 
radio-off phase observations. As a result, we were able to ob- 
tain a distribution of radio-on and -off spin-down rates for each 
data interval, from which we determined estimates for their av- 
erages and associated uncertainties (see Table @. We note that 
the overall average values, that is (u on ) = — 16 ± 1 x 1CP 15 s~ 2 
and (i>off) — —10-8 ± 0- 4 x 10~ 15 s™ 2 , are found to be consistent 
with lKramer et alj <l2006t). 

Following lKramer et al. I ( 120061) , we assume that the spin-down 
in the radio-off states is governed by pure magnetic dipole ra- 
diation. Whereas, we assume that the spin-down in the radio-on 
states is supplemented by an additional torque due to an outflow- 
ing plasma (i.e. a particle wind) and its associated return current. 
Using this model, we can estimate the plasma charge density at- 
tributed to this wind component, to an order of magnitude, via (see 



iKramer et aThood and references therein, for more details): 
3/ Az> 



Ppla 



(3) 



where / is the moment of inertia of the object (taken to be 
10 45 g cm 2 ), Av is the change in spin-down rate of the pulsar be- 
tween emission states (v on — v s), Rpc — \/2irR s v/c is the polar 
cap radius of the pulsar (where the pulsar r adius R is taken to be 
10 6 cm; see, e.g.. lLorimer & Kra mer 2005) and the surface mag- 
netic fiel d strength (for the simplest case of an orthogonal rotator; 
see, e.g., |jackson|[l962l) is 



3.2 x 10 1 



-Vott 

Is- 2 



Gauss . 



(4) 



By substituting R pc and B s into Eqn.[^j|, p p i asma was directly cal- 
culated for each data interval via: 



Ppla 



~ 7.1 xlO 5 



Ai> Vis -1 



ls- 



V2/ ls -2\l/2 



foff 



Cm" 3 . (5) 



These quantities were also compared with the Goldreich- 
Julian charge density of the o bject, which was estimated , to an or- 
der of magnitude, by (see, e.g.. lLorimer & Kramej|2005l ): 



Using the average values for i> g = —10.8 ± 0.4 x 10 -15 s~ 2 and 
v = 1.228965 ± 0.000001 Hz, obtained from the fitting process, 
we estimate p GI to be approximately 0.033 C m -3 . Table|4]shows 
the results of this analysis. We note that different values are mea- 
sured for !> on and i> ff over neighbouring data intervals. However, 
given the uncertainties in these parameters, we consider the differ- 
ent spin-down rates for each emission phase to be consistent with 
the overall averages of i> on and v ff. This is supported by the re- 
sults of an Anderson-Darling test, from which we conclude that the 
distributions of i> on and v a s follow Gaussian statistics. As a result, 
we do not consider the variations in p p i asma with respect to p GJ to 
be significant. 

4.2 Measurement of the braking index 

In the radio-off phase, the rotational slow-down of the pulsar i> e is 
thought to be indicative of the torque produced by magneto-dipole 
radiation dKramer et alj|2003) . Consequently, the rate of change 
in the spin-down rate in the radio-off phase i> g should provide 
a measure of the braking in dex associated with this radiation (e.g. 
iTavlor & Manchested[l977l) 



z> 2 



(7) 



which is not contaminated by the additional electromagnetic torque 
of the plasma current flow; that is, assuming that the additional 
torque completely disappears in the radio-off phase. This is partic- 
ularly important as measurement of this parameter can offer sig- 
nificant information about how the pulsar undergoes its energy 
loss and, ulti mately, provide us with insight into t he electrodynam- 
Xu&Oiaoll200ll : IWu et al.ll2003l : lRudermanll200l: 



original set of fit parameters, to obtain a phase-coherent timing solution due 
to uncertainty in the radio emission phase transition epochs. 



8 Note that an additional factor of 10 7 /c is used to convert between CGS 
units (statC cm -3 ) and SI units (C m -3 ), so as to be consistent with the 



literature (see also Eqn.|6). 
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Table 4. Summary of the results from the residual fitting analysis. The epoch of the spin parameters for each fit is denoted by Pcpoch, T is the length of 
each fit interval and u on , i> Q ff are the radio-on and radio-off spin-down rates respectively. The change in spin-down rate (i.e. i> on — i> a s) is represented by 
Ai>. Definitions for p G1 and p p i asma can be found in the text. Standard 1-cr errors are provided in the parentheses, after the parameters (where possible), and 
represent the least significant digit. 



Pepoch 


T(d) 


u on (10- 15 s" 2 ) 


u oS (10- 15 s- 2 ) 


Au/u^ (%) 


Pplasma (10 2 Cm 3 ) 


Pplasma/Pcj (^°) 


51869.7 


107.5 


-15.2(2) 


-11.13 (9) 


36.6 (6) 


2.5 


73 


52842.3* 


157.5 


-16.0 (2) 


-10.78 (7) 


48.4 (7) 


3.2 


96 


53782.6 


97.7 


-17.6(5) 


-10.4 (2) 


69 (2) 


4.5 


134 


53987.1 


193.4 


-19(1) 


-10.2 (2) 


83 (5) 


5.4 


160 


54340.7 


326.8 


-16.2 (4) 


-10.7 (1) 


51 (1) 


3.4 


102 


54676.6 


133.0 


-17.1 (8) 


-10.3 (3) 


66(4) 


4.3 


130 


54875.0 


196.0 


-15.8(2) 


-11.07 (8) 


42.7 (6) 


2.9 


87 


55198.2 


262.2 


-14.3(1) 


-11.40 (3) 


25.4 (2) 


1.7 


52 



* Pcpoch of data used to derive spin-parameters in lKramer el"ai]|200d ). 



52000 32500 53000 53500 54000 54500 55000 

MJD (day) 

Figure 10. Linear regression fits of o ff (top trace) and u on (bottom trace) 
for PSR B1931+24. Note the anti-correlation in the data, which affects the 
quality of the fits (see text for details). 



down rate in PSR B 1 93 1 +24, by fitting changes in rotational phase 
to the times-of-arrival (using the first three terms of Eqn. 7 in 
Weltevrede et al ■1201 lh . These changes in rotational phase are anal- 
ogous to glitches, that is they are discrete changes in spin-down 
rate, although there are no jumps in spin frequency. In our model, 
each emission phase has a unique spin-down rate (i.e. there is only 
one value for i/ on and z> fj), which we assume is acquired directly 
after (or before) the last (or first) radio-on observation of each ac- 
tive emission phase. Through fitting the observed times-of-arrival 
for changes in spin-down rate only, we find that this method pro- 
duces results consistent with Table|4] However, when the emission 
phase transition times are also taken as a fit parameter we obtain 
very poor fits. We applied maximum and minimum time constraints 
to the 'glitch epochs' but, due to the absence of data bounding a 
given transition into (or out of) a radio-on phase, the fitting process 
could not converge; a global solution could not be achieved because 
times-of-arrival are limited to the radio-on phase data. 



Table 5. The results of the linear regression analysis of the spin-down rate 
data, with least significant uncertainties in the parentheses. The length of 
each fitted data-set is denoted by T. The correlation coefficient and two- 
sided probability value for each fit are given by R and p respectively. The 
rate of change in the spin-down rate is represented by i> of j , with the asso- 
ciated 1-cr error quoted in the parentheses (in units of the least significant 
digit). 



T(d) 


R 


P 


i> ofl (10- 19 s" 3 ) 


3328.5 


0.01 


0.97 


0.05 (156) 


2355.9 


-0.42 


0.35 


-2 (2) 


1415.5 


-0.78 


0.06 


-7 (3) 



Here, we fitted the derived values for i> on and i> ff for a num- 
ber of epochs to determine if significant second order variations 
could be measured (Fig. 1 10b. Unfortunately, however, we could not 
obtain a significant value for p s using data even with the greatest 
observation density (Table|5] see also *j5]for a discussion on this). 
As a result, we could not obtain an accurate measure of the braking 
index of the source. 

Consequently, we sought to determine a value for braking in- 
de x using another method . Using the timing software developed 
bv lWeltevrede et al] J201 ll) . we also modelled the variation in spin- 



4.3 Long-term evolution in spin-down rate 



Using the results from the residual fitting in §4T) we were unable 
to determine a significant value for D s or n for the pulsar. This 
is most likely due to the transition time errors, between emission 
phases, which introduce an anti-correlation between v ori and i> ff (a 
linear fit to these data obtains a correlation coefficient R — —0.9 
and two-sided p-value p — 0.0003; see also Fig. 1 10b. That is, sys- 
tematic over- and under-estimation of the transition times, which 
are not fitted for, cause z> on and z> ff to become anti-correlated 
which, in turn, increases the noise contribution to i> jf . 

To complement this analysis, therefore, we decided to exam- 
ine the long-term evolution in the spin-down rates over the entire 
~ 13-yr data-set. For this analysis, we used timing measurements 
of PSR B 193 1+24 to estimate the contribution of v ori to the rota- 
tional frequency v over neighbouring radio-on phases. We assumed 
that the change in rotational frequency Av, resulting from variation 
between v on and i> ft , can be obtained from fitting the residuals of 
three successive radio-on phases in separate pairs (see Fig. II lb. In 
order to separate the effects of the different spin-down rates in these 
fits, we first form residuals by using a fixed u, set to be the aver- 
age z> ff (using the average determined from residual fitting, i.e. 



(z>off) 



-10.8 ± 0.4 x 10" ia s - ^). We fitted the timing residu- 



als of the latter and first halves of a given pair of radio-on phases, 
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Figure 11. Timing residuals of PSRB 193 1+24, over approximately 80 d of 
data, depicting the Au fitting process. The epochs where u was calculated 
are represented by Pepoch 12 and Pepoch 23 for the first and second pair 
of radio-on phases respectively. The epoch at which Au, hence Au, was 
measured is denoted by Pepoch flt . The bounding regions of the fits are 
shown by the 'start' and 'finish' lines. 



using PSRTIME0, to estimate a value for the rotational frequency 
which is governed by the radio-off spin-down rate. We assumed 
that the contribution of u on to these residuals was negligible due to 
the short length of time the radio-on phases cover (w.r.t. the radio- 
off phases). By obtaining values for u, for the first and second pair 
of radio-on phases, we were able to estimate the total change in ro- 
tational frequency due to the difference in spin-down rate Av over 
the central radio-on phase. Given that u on should only affect the 
residuals over this period, we obtain the relation 



Au 



V23 — V12 



AV X t on : 



(8) 



where Ai> represents u on — u g and t on is the total time spent in the 
central radio-on phase. The spin-frequencies for the first and sec- 
ond fit intervals are U12 and ^23 respectively, as shown in Fig. II II 

To increase the number of measurements, we a nalysed the 
data u sing a stride-fitting method similar to that used bv lLvne et al.l 
we analysed groups of radio-on phase data in steps of sin- 
gle emission phases. We only included data windows which had 
three well defined radio-on phases i.e. containing gaps no greater 
than 10 d between observation^ 1 *! with uncertainties in the emis- 
sion phase transitions less than 4 d (less than 2 d for the longest 
data-set; see Table [6]l. We corrected for the long-term spin-down 
behaviour of the source by subtracting the av erag e spin-down rate 
from the frequency data using the vfit program 11 !. 

Noting that Eqn.[8]is analogous to the equation for a straight- 
line, with a gradient equal to Au and zero intercept, we were able 
to fit Au vs t on across numerous data intervals, as shown in Fig. 1121 
and Table [6] We find that these data are highly linearly correlated, 
which suggests that u on and v a g are well defined. We find that there 
is no evidence for variation in the spin-down rates associated with 
the different phases of emission within our measurement sensitiv- 
ity. 



9 http://www.jb.man.ac.uk/pulsar/observing/progs/psrtime.html 

10 The minimum total time spent in a consecutive radio-on and -off phase 
is 10 d. A gap between residuals of this timescale would cause a large uncer- 
tainty in the contribution of u to Au, due to possible unmodelled rotational 
behaviour. 

11 http://www.jb.man.ac.uk/pulsar/observing/progs/vfit.html 




Figure 12. A linear fit to the change in rotational frequency ( Au) against 
radio-on phase length (t on ) for approximately 13 yr of observations of 
PSRB1931+24. 



Table 6. The linear fit parameters for the long-term spin-down rate study 
of PSR B 193 1+24. The MJD mid-point of each data-set used is denoted 
by MJD m i<j. The total length of each data-set is represented by T and 
N is the number of data points used in each fit. The change in spin-down 
rate obtained from the fitting procedure is denoted by Au and u oll is the 
corresponding radio-on spin-down rate. Standard (l-<x) errors for Au and 
!> on are provided in parentheses after the parameters, and represent the least 
significant digit. 



MJD mid 


T(d) 


N 


Ai>(10- 15 s" 2 ) 


!>on (10- 15 S" 2 ) 


51646.0 


1426 


8 


-5.6(8) 


-16(1) 


53143.4 


1451 


9 


-4(1) 


-16(1) 


54602.8 


1467 


29 


-5.2(3) 


-16.0 (5) 


53356.2 


4688 


45 


-5.2 (2) 


-16.0 (4) 



5 DISCUSSION 

We have shown that PSR B 193 1+24 exhibits regulated emission 
modulation; that is, to the limit of our measurements it maintains 
a constant ADC. To improve on these results, even more frequent 
sampling would be required. This would only be feasible with a 
large dedicated telescope, or new facilities such as LOFAR and the 
sk43 which will have multi-beam capabilities and, hence, greater 
ability to dedicate observing time to individual sources. If we were 
to observe this source over hourly timescales, rather than daily, we 
would dramatically increase our chances of observing the pulsar 
transition between emission modes. As a result, we would possibly 
be able to uncover correlated changes in emission following (or 
preceding) radio-off phases which might give us insight into what 
causes the magnetospheric configuration to change so dramatically. 

We also characterised the modulation timescales of the radio 
emission in PSR B 193 1 +24. We find that the source exhibits a peri- 
odic modulation timescale of approximately 38 d on average. There 
do appear to be variations around this basic periodicity, but it re- 
mains remarkably stable over many years. This degree of stability 
provides a challenge for models of this process as it is significantly 
lon ger than the expected dynamic and plasma timescales (~ ms; 
e.g. lRuderman & Sutherland! 19751) . 



12 Phase 1 science operations with the SKA are projected to commence 
in 2020, with full operation of the telescope proposed in 2024. See also 
http://www.skatelescope.org/about/project/ for more details. 
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Using a residual fitting method, we find that the variations in 
!> on and i> ff , throughout the ~ 13-yr data set, are consistent with 
the measurement uncertainties (see ^ 14. 1 b - This result is supported 
by the long-term analysis (see § [43}, which provides additional 
evidence to suggest that we do not observe any significant changes 
in the spin-down rates attributed to each emission mode; that is, 
the pulsar appears to retain a constant Ai> = 48 ± 2% between 
emission phases. This implies that the apparent changes in plasma 
flow in the radio-on phase are also not significant ((/9 p i a sma) = 
4 ± 1 x 10" 2 Cm -3 ; see Table |4](, thus indicating a surprisingly 
high degree of stability in the bi-modal system. 

We note, however, that there is an unavoidable systematic ef- 
fect in the short-timescale fitting procedure. More specifically, the 
errors in the transition times between emission phases introduced 
uncertainty into the calculation of i/ on and i> g. This is highlighted 
by the results in Fig. QJJ] where the interdependency between the 
two spin parameters is made evident. We see that these quanti- 
ties are in fact almost perfectly anti-correlated. In our model, the 
transition times between phases of emission are not taken as a fit 
parameter. Therefore, observations spaced farther apart between 
two emission modes will have a higher uncertainty in the transi- 
tion time compared with two spaced closer together. Furthermore, 
if the amount of time spent in a radio-on phase is overestimated, 
then the amount of time in a consecutive radio-off phase will be 
underestimated accordingly. Given that radio-on phases are much 
shorter than radio-off phases on average, it is evident that the effect 
on the error of v orL will be greater than that of z> D ff . This naturally 
explains the observed anti-correlation between the two parameters 
and the larger error bars in i/ on . Therefore, it is possible that we 
did not see any significant variation in the spin-down rates due to 
this systematic effect. To overcome this problem, we again would 
require greater observation cadence to provide better constraints on 
the transition times. 

We also attempted to obtain a value for the braking index of 
PSR B 193 1+24 from v a, in the hope of determining information 
about the energy loss of the object. Using Eqn.[7] we predicted the 
i) that would be required to obtain a braking index n = 3. As- 
suming the average parameters for v and v a g from residual fitting, 
we would expect v = 2.85 ± 0.05 x 10~ 28 s " 3 (i.e. ~ 10 9 less 
than the upper limit derived from fitting separate residual data-sets). 
Given the level of timing noise in this pulsar, and the v quoted in 
iHobbs et al.l ^OlOlP^l it seems very improbable that any increase 
in data density would allow us to obtain a value for v as low as that 
predicted from theory. The problem lies in decoupling the long- 
term contribution of z> on from the timing residuals; that is, the value 
for v from global fitting will always be contaminated by the modu- 
lation of v. In the case of local fitting (!> on , off), there will probably 
always be an uncertainty in the transition times which, again, means 
that it will be difficult to decouple z> on , off from !> av . 

With the next generation of telescopes (e.g. LOFAR, FAST, 
ASKAP, MeerKAT and the SKA) coming online in the near fu- 
turj^l we anticipate a dramatic increase in the number of known 
intermittent pulsar f^l Coupled with the continued observation of 
known intermittent sources (e.g. PSRs B0823+26, J1832+0029, 

13 IHobbs et al l feOlOl) obtain v = 7.0±0.1x lCT 25 s " 3 for a time span 
of ~ 13.7 yr, which is a factor of 10 6 less than the upper limit derived from 
this work. As they apply a global fit to the timing data, this lower limit is 
corrupted by the spin-down rate in the radio-on phase. 

14 Although LOFAR is already collecting data, the complete array is yet to 
be fully built. 

15 Approximately 900 pulsars are predicted to be discovered with LOFAR 



J1841-0500 and B1931+24), characterisation of these objects 
should lead to an improved understanding of the general emission 
and rotational properties of pulsars which, ultimately, should facil- 
itate the development of more realistic radio emission and magne- 
tospheric models. These observations should also provide further 
insight into how the different 'types' of nulling pulsars, i.e. normal 
nulling pulsars, RRATs and longer-term intermittent pulsars, are 
related. 



6 CONCLUSIONS 

In this work, we have confirmed an d expanded substantially upon 
the findings of iRramer et alj d2006l) . which has enabled us to pro- 
vide furth er evidence for the magnetospheric - state switching sce- 
nario (e.g. lBartel et al.lll982l ; Evne et al]|2010l ; lTimokhinll2010|) . In 
order to determi ne how these events are tri ggered (e.g. via 'circum- 
pulsar asteroids' ; ICordes & Shannonl2008l) . however, it is clear that 
dedicated infra-red and high-energy observations of this source are 
required. 

If the spin-down rates attributed to each emission mode, and 
ADC, for this pulsar are truly constant over time, as suggested 
by our observations, this would be quite remarkable; as of yet, 
there is no clear reason to suggest why a pulsar should retain a 
memory of its previous magnetospheric state i.e. particle flow (c.f. 
iLi et alj[2012l) . This raises a couple of important questions; 1) Is 
there charge or matter transfer in the magnetosphere which leads 
to these separate regulated states and, if so, how would this occur? 
2) Why does the pulsar consistently assume the same rotational 
and emission characteristics? These questions are likely only to be 
answered based on future observations of this source and others 
like it, which may as yet be undiscovered. We note, however, that 
the rate at which new intermittent sources sources are discovered 
should increase significantly with the use of next generation tele- 
scopes. Therefore, we may only need to wait a few years before 
significant breakthrough is made in magnetospheric modelling and 
our understanding of emission cessation in pulsars. 
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